Gravitational Theory, Galaxy Rotation Curves and 
Cosmology without Dark Matter 



J. W. Moffat 

The Perimeter Institute for Theoretical Physics, Waterloo, Ontario, N2J 2W9, 

Canada 

and 

Department of Physics, University of Waterloo, Waterloo, Ontario N2Y 2L5, 

Canada 

Abstract 

Einstein gravity coupled to a massive skew symmetric field F^ux leads to 
an acceleration law that modifies the Newtonian law of attraction between 
particles. We use a framework of non-perturbative renormalization group 
equations as well as observational input to characterize special renormaliza- 
tion group trajectories to allow for the running of the effective gravitational 
coupling G and the coupling of the skew field to matter. Strong renormal- 
ization effects occur at large and small momentum scales. The latter lead to 
an increase of Newton's constant at large galactic and cosmological distances. 
For weak fields a fit to the flat rotation curves of galaxies is obtained in 
terms of the mass (mass-to-light ratio M/L) of galaxies. The fits assume that 
the galaxies are not dominated by exotic dark matter and that the effective 
gravitational constant G runs with distance scale. The equations of motion 
for test particles yield predictions for the solar system and the binary pulsar 
PSR 1913+16 that agree with the observations. The gravitational lensing of 
clusters of galaxies can be explained without exotic dark matter. An FLRW 
cosmological model with an effective G = G{t) running with time can lead to 
consistent fits to cosmological data without assuming the existence of exotic 
cold dark matter. 

e-mail: jmoffat@perimeterinstitute.ca 

1 Introduction 

A nonsymmetric gravity theory (NGT) has been apphed to explain rotation curves 
of galaxies and cosmology, without invoking dominant dark matter and identifying 
dark energy with the cosmological constant [HIS El El- In the following, we develop 
a simpler gravitational theory based on Einstein's general relativity (GR) and a 
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skew symmetric rank three tensor field -F^i/a, forming a metric-skew-tensor-gravity 
(MSTG) tlieory. We sliall apply tliis gravity tlieory to explain the fiat rotation 
curves of galaxies and cluster lensing without postulating exotic dark matter. The 
weak field approximation of the NGT field equations corresponds to the equivalent 
field equations in MSTG and, therefore, to this order of approximation will yield the 
same predictions as MSTG. The more complicated NGT will have different possible 
consequences for strong gravitational fields. However, in the following, we shall be 
concerned with the observational data for the solar system, the binary pulsar PSR 
1913+16, galaxies, clusters of galaxies and cosmology. The MSTG theory provides 
a more direct way of deriving a comparison with the data. 

Since no dark matter has been detected so far, it seems imperative to seek a 
possible modified gravitational theory that could explain the now large amount 
of data on galaxy rotation curves. A cosmological model obtained from the field 
equations and a running of the effective gravitational coupling constant G in the 
matter epoch could explain the growth of large scale structure formation without 
invoking cold dark matter. The running of the cosmological constant would pro- 
duce a quintessence-like dark energy that could account for the acceleration of the 
expansion of the universe |Sl El E! • 

A renormalization group (RG) framework [H] for MSTG is developed to describe 
the running of the effective gravitational coupling constant G, and the effective cou- 
pling constant 7c that measures the strength of the coupling of the F^^x field to 
matter. A momentum cutoff identification k = k{x) associates the RG scales to 
points in spacetime. For the weak field static, spherically symmetric solution the 
RG flow equations allow a running with momentum k and distance i = 1/k for the 
effective Newtonian coupling constant G = G{r), the coupling constant 7c = 'jdr), 
and the effective mass of the skew field fi = fi{r) where r denotes the radial coordi- 
nate. The form of G{r) as a function of r, obtained from the modified Newtonian 
acceleration law, leads to agreement with solar system observations, terrestrial grav- 
itational experiments and the binary pulsar PSR 1913-1-16 observations, while strong 
renormalization effects in the infrared regime at large distances lead to fits to galaxy 
rotation curves. 

A fit to both low surface brightness and high surface brightness galaxies is 
achieved in terms of the total galaxy mass M (or M/ L) without exotic dark matter. 
A satisfactory fit is achieved to the rotational velocity data generic to the elliptical 
galaxy NGC 3379. Fits to the data of the two spheroidal dwarf galaxies Fornax and 
Draco and the globular cluster u Centauri are also obtained. The predicted light 
bending and lensing can lead to agreement with galaxy cluster lensing observations. 
A model of the modified acceleration law that includes a description of radial ve- 
locity curves in the core of galaxies as well as in the outer regions of the galaxy is 
shown to yield good fits to rotational velocity data. 

For a homogeneous and isotropic universe, based on the Friedmann- Lemaitre- 
Robertson- Walker (FLRW) model, the average value of the skew symmetric field 
Ffj^i^x is zero. However, for an effective gravitational constant G = G(t) that runs 
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with time, we find that the constraint G{t) ~ Go (where Go is Newton's constant, 
Go = 6.673 X 10^^ g^^ cm^sec"^) must be imposed on G{t) at the time of big bang 
nucleosynthesis (BBN) and the surface of last scattering (decoupling) when the 
fraction of baryons is fi^ ~ 0.04. For a rapid increase of G{t) after the time of 
decoupling when G{t) reaches a value G{t) ~ 6Go, the baryonic matter fraction 
becomes Qm ~ G^b ~ 0.24 and there is no need for undetected non-baryonic dark 
matter^ 



2 The Field Equations 

The action for the MSTG theory is given by 

S = Sg + Sf + Sfm + Sm, (1) 
where Sg is the Einstein-Hilbert action (we choose units such that c = 1): 

Sg = / d'x^iR - 2A), (2) 

where g = Det{gfj,u), g^y is the symmetric metric tensor of pseudo-Riemannian geom- 
etry, R = g^^R^u is the Ricci scalar and A is the cosmological constant. Moreover, 
Sp is the skew field action 

Sp = j d'xV^{^F^y,F^''P - ^f^'A^^An, (3) 

where 

d^^A^x + d^Axf, + dxA^y, (4) 

and fi is the mass of the skew field. The gravitational constant G in the action Sg 
is defined in terms of the "bare" gravitational constant Gq: 

G = GoZ, (5) 

where Z corresponds to a "renormalization" of G. 
The actions Sm and Sp satisfy the relations 

1 6Sm 1 1 6Sf 1„ 



2 ' V^Sgt"' 2 
Moreover, we have 

= _j (7) 

^An alternative cosmological model has been proposed [H| that fits the data without a cosmo- 
logical constant and for a universe consisting only of baryons and cold dark matter and with a low 
value of the Hubble constant: Ho ^ 46 km/s/Mpc. 
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Here, TMfn^ is the energy-momentum tensor for matter, while Tpfj.^, is the energy- 
momentum tensor containing the contributions from the massive skew field A^i^. 
Also, Jf^i, is the tensor density source for the A^^ field. 

A possible action for the skew field matter coupling is pijj: 

SpM = I rf"xA„.J*^^^ = -3 / rf"xe"'^'^M.«a, J, 



FM = J U -I^^Xixv-J = -O J U Xt ■ /1q,/3C^J,,, [O) 

where is a vector current that can be conserved and J*'^'^^ = e^'^^^J^ is the dual 
tensor current density. The Levi-Civita symbol e^"^'^ transforms as a contravariant 
tensor density, while e'^'^^°' = ±1, 0. The Pauli coupling (jSJ is gauge invariant for any 
J^. The source current can be identified with a matter fermion current associated 
with baryon and lepton number. 

The field equations derived from the action principle are given by 

G^, + Ag^, = SirGT^,, (9) 
V^F^.. + fx^A^, = -L (10) 



Ill/a I -''■III/ ; '-'I1U1 

where Gf^u = R^iu — ^QfiuR, T^u = Tm^u + Tp^u and V°" denotes the covariant 
derivative with respect to g^^. Moreover, we have 

J^i. = e^^^pcTJ^. (11) 

The Bianchi identities satisfied by the Einstein tensor G^'^ lead to the conserva- 
tion laws 

V^T^'^ = 0. (12) 



3 Equations of Motion of Test Particles 

The equations of motion for a test particle are given by 

^ + |^^^}«"n'^ = ^?^"/..n^ (13) 

where r is the proper time along the path of the particle and = dx^/dr is the 
4-velocity of the particle. Moreover, 

{ /il } " ^^^^'"^ ^ ^^^''^ ~ ^''"'''^ ' ^^^^ 
is the Christoffel connection and 

fail = A(9[q, ( H„yxg^l]r, I • (15) 



4 



Here, H^,yx is given by 

H^ux = ^(5aV + ^M^-A + d^Ax^), (16) 

and A is a coupling constant with the dimension of length that couples the skew 
field to the test particle. 

For a spherically symmetric, static skew symmetric potential field A^^, there are 
two non- vanishing components, the "magnetic" field potential AQr{r) = w{r) and 
the "electric" potential field Aq^^t) = f{r)sm9. We shall assume that only the 
electric field contribution /(r)sin6' is non-zero (no static magnetic poles). Then, 
the tensor -ff^i/A has only one non- vanishing component: 

He^r = ^drAe^ = f sine, (17) 



where /' = df /dr. We obtain 



^ dr \^/crfr^ 



;i8) 



For a static spherically symmetric gravitational field the line element is 

ds"^ = ^{r)dt^ - a{r)dr'^ - r'^{de'^ + sin^ ed<i?). (19) 

The equations of motion for a test particle are given by 

(iV a' fdry r/dOy fsm^ 9\fd(l)\^ 'y' / dty 
dr"^ 2a Kdr J a Kdr J \ a J Kdr J 2a Kdr J 

~^ adr (^a7r^) (rfr) ' ^ ^ 

d^^t(^) + l±(JnL_) = 0, (21) 

dr"^ 'J Kdr J Kdr J 'j dr\y/cryr^J KdrJ ' 

The orbit of the test particle can be shown to lie in a plane and by an appropriate 
choice of axes, we can make 9 = 7i/2. Integrating Eg. (j23|l gives 

r'f = J, (24) 

dr 

where J is the conserved orbital angular momentum. Integration of Eq. ()21|1 gives 

dt 1 



dr 7 



y/ctyr^ 



(25) 
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where E is the constant energy per unit mass. 

By substituting into and using (j211), we obtain 

d'r a' fdry J' j / X^f ^ ^ 1 d / \ f X^f 



We do not have an exact spherically symmetric, static solution for massive 
MSTG. For large enough values of r, we shall approximate the metric components 
a and 7 by the Schwarzschild solution: 

~ ^ 2GAf ^ 7W ~ 1 (27) 

r 

and make the approximations that 2GM/r <^ 1, Xf /r"^ <^ 1, //r^ ^ 1 and the slow 
motion approximation dr/dt -C 1. Then, for material particles we set = 1 and 
becomes 



dt"^ dr V 



where Jn is the Newtonian orbital angular momentum. 

4 Linear Weak Field Approximation 

We expand gfj_u about a Minkowski flat spacetime 

g^,^ = r]^^ + hf,^ + 0{h'^), (29) 

where 77^,^ = diag(l, — 1, — 1, — 1) is the Minkowski metric tensor. The skew field 
-^/ii/A obeys the linearized equations of motion for empty space 

d^F^,^ + fx'A^, = 0. (30) 

These are the massive Kalb-Ramond-Proca equations which are free of ghost 
pole instabilities and possess a positive Hamiltonian bounded from below. 

For the components Ag^ = f{r) sin6, we find from (|3U|) for a static spherically 
symmetric field: 

/»--/'(r)-/iV(r)=0, (31) 
r 

which has the solution 

/(r) = -sG^M^ exp{-nr){l + fir), (32) 
3 

where s is a dimensionless constant. We note that Agfj^ are the components of a 
second rank skew tensor, which are not equivalent for a massive skew field to a 
pure scalar field. We obtain (jHT|) when we transform to the spherically symmetric 
equations in polar coordinates. 
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Let us now choose g^^ to be a Ricci-flat GR background metric Qfj^y = g 



GR 
fiu ■ 

The skew field F^^^ obeys the hnearized equation of motion in the GR background 
geometry 

V^F,,,, + AA'^^Bp^,, + ^?A^, = 0, (33) 

where V'^ and Bp^^j^ denote the background GR covariant derivative and curvature 
tensor, respectively. 

The skew field components A^r have the solution [T^ : 

{2a" - /i2)Ao. = 0, (34) 

where 7(r) ~ l/a(r) ~ l — 2GM/r. We conclude from this that Aq^ vanishes outside 
the source. In the massless case, this is the surviving spherically symmetric ghost 
mode which in this case is pure gauge. When the mass term is added, although 
these modes are no longer pure gauge, they do not propagate since they are locally 
coupled to the source. Thus, in the absence of sources only the Ag^ components are 
physically interesting. 

For the static, spherically symmetric spacetime, the linearized equation of motion 
for /(r) on a Schwarzschild background takes the form 

^ 2GM\j^„ 2(^ 3GM\^, f /j,^ , 8GM 



/"-- 1-^ /'-^ + ^ / = 0, (35) 
/ r \ r / \r"^ r / 

The solution to this equation in leading order is [U [T21 HHI : 

fir) = '^3^' T\~gm (l + /ir + ^ [2 + /xr exp(2/.r)Ez(l, 2/.r) (/ir - 1)] ) , (36) 



(/xr) 

where Ei is the exponential integral function 
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Eiin,x)= I dt'-^^^y^. (37) 



We obtain for large r the solution 



l 5G^M^exp(-/ir)(l+/ir) 

/(O = 3 in^^ • ^''^ 



For fiGM <^ 1, we have (fir)'^^^^ ~ 1, yielding the solution 

To the order of weak field approximation, we obtain from Eqs. ()28|) and 



^"^ + "^"^^-^^^ (1 + ^, (39) 



where the constant a is given by 



(40) 
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We have required in Eq. ()39|l that the additional acceleration on the right-hand side is 
a repulsive force. This is in keeping with the identification of A^i, with the potential 
field of a massive axial vector spin 1+ boson. 

An analysis of the physical properties of a massive spin l"*" boson exchange was 
carried out by Damour, Deser and McCarthy (DDM) JU] , based on a modification of 
a massive NGT with skew field B^^^. In a version of a rigorous massive NGT jHlll] , the 
weak field solution for a massive Kalb-Ramond field ^T] follows consistently from a 
weak field approximation to the rigorous NGT field equations and leads to a stable 
vacuum without ghost energy modes. The solution of the weak field approximation 
NGT field equations yields the solution (jH^ for a spherically symmetric system. 

The macroscopic coupling action for the DDM skew symmetric B^^iy field is given 
in their notation by 

Sj = -\fcj d^xHx.ur^^'' = \h J d^xe^''"'^B^pd,J,, (41) 

where fc denotes a dimensionless coupling constant, H\^y is the fully skew symmetric 
field strength and is a current vector that can be conserved. Moreover, J*^^^'^ = 
^Xfiuaj^ is a totally antisymmetric tensor density. The action Sj (within trivial 
numerical factors) is the same action (jH)) that we have used to obtain our field 
equations (jl(J|) for the A^,^ potential field and the equations of motion of a test 
particle (jl3j) . DDM show that the lowest order matter coupling of the skew B-field 
corresponds to a vector particle "fifth force" coupled to the current J^, which can be 
identified with a fermion current with the dimensionless coupling = y/ATrGofifc- 
Note that the coupling constant g^, is inversely proportional to the range of the fifth 
force. 

We can consider the strength of the skew field coupling by introducing the cou- 
pling constant 

a,= A^2/, xlO-i^cm, (42) 
rriN 

where ttt-at is the nucleon mass (one atomic mass unit). Moreover, Oc has dimensions 
of a length and couples baryon number to mjsiJ^. 

The coupling constant ac can be expressed in terms of a^ = g\l (47rm^) as 



Oic = - — = ^/a5ro, (43) 

where ro = 1/yU is the range of the force. This show that ac can take large macro- 
scopic values provided that the range tq is large enough. In our RG flow framework, 
the coupling constant ac will run and become increasingly large in the IR momen- 
tum limit A; — > 0, corresponding to increasing values of the range vq. We can define 
a "renormalized" coupling constant 

ac = acoA, (44) 
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where Oco is the bare skew field couphng constant and A is a renormahzation con- 
stant. 

Experiments have put stringent bounds on possible fifth forces and the magni- 
tude of (75 Pl]. However, DDM make the important points for phenomenological 
purposes that the strength of the coupling fc is unbounded as the range increases, 
and that the magnitude of the B-field (in our notation the B-field is A^iy) is pri- 
marily proportional to fc and independent of the range tq. This allows the B-field 
to have a "gravitational" strength, while still keeping compatibility with the exist- 
ing stringent bounds on possible violations of the weak equivalence principle and 
composition-dependent effects in Newtonian gravity. 

It should be noted that significant observational bounds on fifth forces only apply 
to distances < 100 AU in the solar system and outer solar system, corresponding to 
fi ~ 10^^° eV. The experiments do not place useful bounds on the strength of the 
coupling of a gravitational fifth force at galactic distance scales ~ 5 — 100 kpc or 
at cosmological scales. This is in keeping with the requirement that our modified 
acceleration law will be approximately Newtonian at distance scales within the solar 
system, and that the MSTG predictions for solar system distance scales agree with 
the solar system observations. 



5 Orbital Equation of Motion 

We set 6* = 7r/2 in (fT^ . divide the resulting expression by dr^ and use Eqs. fl^ and 
(OKI) to obtain 



2 J2 

+ 



'•' +E 



E 

a 



(45) 



We have ds^ = Edr"^, so that ds/dr is a constant. For material particles E > and 
for massless photons E = 0. 

Let us set M = 1/r and by using (j211), we have dr/dr = —Jdu/dcj). Substituting 
this into (j45j) . we obtain 



/ du\' 

km) 



E + 



E 



(46) 



By substituting (j27jl and dr/dcp = —{l/u'^)du/d(j) into (jlHl), we get after some ma- 
nipulation: 



d'^u 



+ u 



EGM EXsG^M^ 



J2 



exp 



rou 



1 + 



^) 



3GMu\ (47) 



where tq = l//i. 

For material test particles E = 1 and we obtain 



d'^u 



u 



GM 



3GMu^ 



4expf — ^Vl + — ^ 



J2 



roU 



TquJ 



(48) 



9 



where K = XsG'^M'^ /Sr^. On the other hand, for massless photons ds"^ = and 
E = and dHj) gives 

f^+u = 3GMu\ (49) 

6 Galaxy Rotational Velocity Curves 

A possible explanation of the galactic rotational velocity curves problem has been 
obtained in NGT jTHI. We shall now obtain an equivalent explanation of the rota- 
tional velocity curves from the present MSTG theory. From the radial acceleration 
derived from ()39|) experienced by a test particle in a static, spherically symmetric 
gravitational field due to a point source, we obtain 

GooM exp(-r/ro)/ r\ 
= ^ + ^ 2 U + — ■ (^0) 

Here, Goo is defined to be the effective gravitational constant at infinity 



where Go is Newton's "bare" gravitational constant. Moreover, M is the total mass 
and Mq is a parameter related to the strength of the coupling of the skew field to 
matter. 

This conforms with our definition of G in Eq.(j3)), which requires that the effec- 
tive G be renormalized in order to guarantee that ()50|) reduces to the Newtonian 
acceleration 

O-Newton — — yp'^) 

at small distances r <C tq. The constant a is given by (we reinstate c): 



The integration constant s in is dimensionless and can be modelled as 

s = gM\ (54) 

where g is a coupling constant and 6 is a dimensionless constant. We choose b = —3/2 
and set XgGo/Sc^rQ = \/Mq. The choice of 6 = —3/2 yields for galaxy dynamics 
an approximate Tully-Fisher law ^Sl- The expression (jH^ for the parameter s is 
a phenomenological model for the coupling of the skew field to matter, which does 
not at present have a fundamental interpretation. 
We obtain the acceleration on a point particle 

a(r) = + G^^ r"^-:''"^ (l + -). (55) 
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By using (j5T| . we can express the modified acceleration in the form 



a(r) 



GqM 



1 + 



'Mo 
M 



exp(-r/ro) ( 1 + — ) 



(56) 



In section 10, we shall explain the running with distance scale of the effective 
gravitational constant G = G{r) and the skew field coupling constant s = s{r) in 
an RG flow, effective MSTG action framework in which RG phase space trajectories 
possess Gaussian and non-Gaussian fixed points. 

We can rewrite (jSUj) in the form 



a(r) 



G{r)M 



(57) 



where 



G(r) = G( 



1 + 

V M 



exp( 



-r/ro) IH 

ro 



(5J 



Thus, G{r) describes the running with distance of the effective gravitational constant 
in the RG flow scenario. Moreover, the parameter Mq, which is a phenomenological 
description of the coupling strength of the skew field to matter will also run with 
distance in the RG flow scenario. 

We apply to explain the flatness of rotation curves of galaxies, as well as 
the approximate TuUy-Fisher law ^HI, oc GqM oc L, where v is the rotational 
velocity of a galaxy, M is the galaxy mass 



M = M^ + Mhi + Mdb + M 



(59) 



and L is the galaxy luminosity. Here, M^:,Mhi, Mdb and Mf denote the visible 
mass, the mass of neutral hydrogen, possible dark baryon mass and gas, and the 
mass from the skew field energy density, respectively. 

The rotational velocity of a star v obtained from v'^{r)/r = a(r) is given by 



iGnM 



'Mo 
M 



1 - exp(-r/ro)( 1 + 



ro 



1/2 



(60) 



Let us define the relation between parameters Mq and ro in terms of the magni- 
tude of the constant acceleration 



ao 



GoMo 



(61) 



We assume that for galaxies and clusters of galaxies this acceleration is determined 
by 

ao = cHo. (62) 
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Here, Hq is the current measured Hubble constant Hq = 100/ikms ^ Mpc ^ where 
h = (0.71 ± 0.07) [IE]. This gives 

ao = 6.90 X 10"*cms~l (63) 



We note that oq = cHq ~ {yA/3)c^, so there is an interesting connection between 
the parameters Mq, tq and the cosmological constant A. 

A good fit to low surface brightness and high surface brightness galaxy data is 
obtained with the parameters 

Mo = 9.60 X 1O"M0, ro = 13.92 kpc = 4.30 x 10^^ cm (64) 

and M (or the mass-to-light ratio M/L). In (jMj) the parameter tq is not independent 
of Mq. By using ^ and substituting Mq = 9.60 x IO^Mq into dHD), we obtain 
the ro in (I64|l . Thus, we fit the galaxy rotation curve data with one parameter Mo 
and the total galaxy mass M. Since we are using an equation of motion for point 
particle sources, we are unable to fit the cores of galaxies without supplementing 
the acceleration formula with a galaxy core model based on a mass distribution. 

Let us now describe a model of a spherically symmetric galaxy with a core density 
of visible matter, Pc{r), within a core radius r < r^. The acceleration law takes the 
form 



«W = -2 U + 




l-exp(-r/ro)(^l + ^)]|, (65) 



where ^ 



M{r) = Atc / c/rV>e(r') (66) 





is the ordinary matter inside the luminous core of the galaxy described by a ball 
of radius r = r^. Inside the core radius the dynamics is described by Newtonian 
theory. Moreover, the acceleration outside the core radius r > Tc is described by 
fj56|) . while M is determined by (jSHl)- A simple model for Ai{r) is given by^ 

M{r) = Mi^-^y, (67) 

where M is the constant total mass and (3 = constant. The rotational velocity 
derived from the acceleration law is 



V 




exp(-r/ro) ( 1 + — ) 



1/2 

(6J 



^We can choose to solve a Poisson equation for the potential that gives a more accurate descrip- 
tion of the core behavior of the galaxy, including bulge effects and other detailed features of the 
galaxy associated with a measured luminosity distribution. However, the model considered here 
yields a reasonable description of the velocity curves in the core regions of the galaxies. 
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The modified acceleration law (j65|) can be compared to the Newtonian law using 



The fits to the galaxy rotation curves v in km/s versus the galaxy radius r in 
kpc are shown in Fig. 1. The acceleration law is given for point particles by (j56|l . 
The data are obtained from ref. [TrJ. 

In Fig. 2, fits to two dwarf galaxies (dSph) are shown. We assume that the 
relation between the velocity dispersion a and the rotational velocity v takes the 
simple form in e.g. an isothermal sphere model for which v ~ y2cT. The error bars 
on the data ^Hl for the velocity dispersions are large, and in the case of Draco, due 
to the small radial range 0.1 kpc < r < 0.6 kpc, the Newtonian curve for 



cannot be distinguished within the errors from the MSTG prediction. However, it 
is noted that the MSTG prediction for v appears to flatten out as r increases. For 
Draco M/L = 28.93 + 5O.3O(9.58)(M0/L0), whereas for Fornax M/L = 1.79 + 
O.72(O.4O)(M0/L0). There is also an expected large error in the distance estimates 
to the dSph. Another serious potential source of error is that it is assumed that 
dSph galaxies are in dynamical equilibrium. The two studied here are members of 
the Local Group and exist in the gravitational field of a larger galaxy, the Milky 
Way. Thus, the tidal interactions with the larger galaxy are expected to affect the 
dynamics of dSph galaxies and the interpretations of velocity dispersions [201 ■ These 
issues and others for dSph galaxies are critically considered in the context of dark 
matter models by Kormendy and Freeman i2l^. 

We have also included a fit to the elliptical galaxy NGC 3379 as shown in Fig. 3. 
The elliptical galaxy NGC 3379 has been the source of controversy recently |22]- The 
velocities of elliptical galaxies are randomly distributed in the galaxy. However, the 
gravitational potential that would be experienced by a test particle star or planetary 
nebula in circular rotation about the center of the galaxy can be extracted from the 
line-of-sight velocity dispersion profiles. The data for R/Res > 0.5 refer to planetary 
nebula. 

We use the mean values of the extracted rotational velocities for NGC 3379, 
obtained by Romanowsky et al. [2'2\ and find that our predicted rotational velocities 
agree well with their data as shown in Fig. 3. According to Romanowsky et al. 
there appears to be a dearth of dark matter in the elliptical galaxy which needs 
to be explained by dark matter models and N-body cosmological simulations. For 
Milgrom's MOND [231121123 it is argued by Sanders and Milgrom [23 that NGC 
3379 is marginally outside the MOND regime with an acceleration a ~ (cio)Miigrom ~ 
1.2 X 10~^ cms~^, so the MOND prediction for elliptical galaxies should be closer to 
the Newtonian-Kepler prediction for the rotational velocity curve. Romanowsky et 
al. also give data for the two elliptical galaxies NGC 821 and NGC 4494, but the 



(iniD: 



C'-Newton V 



GoMjr) 



(69) 




(70) 
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intrinsic circular velocities inferred from the line-of-sight velocity dispersion profile 
data are not given by the authors, although the trends of the data are similar to 
NGC 3379. 

A fit to the data for the globular cluster uj Centauri is shown in Fig. 3. The data 
is from McLaughlin and Meylan j|27^. We use the velocity dispersion data and assume 
that the data is close to the rotational velocity curves associated with the velocity 
dispersion dp i.e. the isothermal sphere model relation v ~ \/2ap holds. The fit to 
the data reveals that the predicted rotational velocity cannot be distinguished from 
the Newtonian-Kepler circular velocity curve within the orbital radius of the data. 
The authors conclude that there appears to be no room for dark matter, whereas 
the our results agree well with the data. For Milgrom's MOND, the magnitude of 
the acceleration is larger than the MOND upper limit (ao) ~ lO^^cm /sec and we 
expect to obtain a Newtonian-Kepler rotation curve. 

In Fig. 4, we show fits to the four galaxies NGC 1560, NGC 2903, NGC 4565 
and NGC 5055 using the modified acceleration law including the galaxy core 
with (3 = 2. As can be seen from the data fitting, the agreement with the core and 
extended rotation curves data |2H1 1211 is good for the two fitting parameters M and 
the core radius Tc- 

In Fig. 5, we display a 3-dimensional plot of v versus the range of distance 
0.1 kpc < r < lOkpc and the range of total galaxy mass M used in the fitting 
of rotational velocity data. The red surface shows the Newtonian values of the 
rotational velocity f , while the dark surface displays the prediction for v obtained 
from our gravity theory. 

Table 1 displays the values of the total mass M used to fit the galaxies and the 
mass-to-light ratios M/ L estimated from the data in references given in J7j , using 
the non-core acceleration model ()56p . 

In Milgrom's phenomenological model j23| |23 we have 



where (ao)Miigrom = 1-2 x 10~'^cms~^. We see that (f7T|) predicts that the rotational 
velocity is constant out to an infinite range and the rotational velocity does not 
depend on a distance scale, but on the magnitude of the acceleration (ao)Miigrom^- 
In contrast, our modified acceleration formula does depend on the radius r and 
the distance scale tq which for galaxies is fixed by the formula ()62j). The MSTG 
velocity curve asymptotically becomes the same as the Newtonian-Kepler prediction 
as r ^ cxd: 



where Goo is the renormalized value of Newton's constant. 

Using the Sloan Digital Sky Survey (SDSS), Prada et al. [301 have studied the 
velocities of satellites orbiting isolated galaxies. They detected approximately 3000 

^In formulations of Milgrom's model, the modified Newtonian potential has a logarithmic de- 
pendence on r 



v'^ = GoM(ao)Mii! 



Igrom; 



(71) 




(72) 
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satellites, and they found that the line-of-sight velocity dispersion of satellites de- 
clines with distance to the primary. The velocity was observed to decline to a 
distance of ~ 350 kpc for the available data. This result contradicts the constant 
velocity prediction (HH) of MOND, but agrees with the MSTG prediction (f72|) . It also 
agrees with the cosmological models which predict mass profiles of dark matter halos 
at large distances. During the last two decades of numerical modelling of galaxy 
formation, they have produced a density profile of dark matter halos, p oc at 
large radii, which does not depend on the nature of the dark matter pilj . 

Although it appears that the MOND prediction for the velocity distribution of 
satellite galaxies disagrees with the Prada et al. data, it can be argued that if a 
system is embedded in an external gravitational field of Newtonian acceleration Oext, 
then if the internal Newtonian acceleration of the system aNewton satisfies ONewton < 
Oext < flo, then the modified Milgrom acceleration is aNewton (ao/^ext) instead of 
i/aociNewton- The acceleration becomes Kepler-Newtonian and anisotropic with an 
increased effective gravitational constant |26j. However, the MSTG prediction of 
the velocity dispersion does not depend sensitively on knowing the modelling of an 
external gravitational field influence on the satellite galaxies. 



7 Local and Solar System Observations 

We obtain from Eq. ()48p the orbit equation (we reinsert the speed of light c): 



(fu GM K , , , 



1 + 



+ — ^« , (73) 



c2 



where now K = XsG'^M'^ /Sc^Tq. Using the large r weak field approximation, and 
the expansion 

exp(-r/ro) = l-^ + ^(^) + ... (74) 
we obtain the orbit equation for r ^ tq: 

+ u = A^ + 3^M^ 75 

where 

^=^-12- (76) 

We can solve Eq. ()75|) by perturbation theory and find for the perihelion advance 
of a planetary orbit 

A^= (GMe-c^i^e), (77) 

C Li 

where Jat = {GMqL/ c^Y^"^ ^ L = a(l — e^) and a and e denote the semimajor axis 
and the eccentricity of the planetary orbit, respectively. 
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We now use the running of the effective gravitational couphng constant G = 
G{r), determined by (jKH|l and find that for the solar system r -C tq = 14kpc, 
we have G ~ Gq within the experimental errors for the measurement of Newton's 
constant Gq. We choose for the solar system Kq -C 1.5 km and use G = Gq to 
obtain from ()77|) a perihelion advance of Mercury in agreement with GR. 

For terrestrial experiments and orbits of satellites, we have also that G ~ Go for 
r ~ where t"® <C = 14 kpc is the radius of Earth. We then achieve agreement 
with all gravitational terrestrial experiments including Eotvos free-fall experiments 
and "fifth force" experiments. 

For the binary pulsar PSR 1913+16 the formula (f77j) can be adapted to the 
periastron shift of a binary system. Combining this with the MSTG gravitational 
wave radiation formula, which will approximate closely the GR formula, we can 
obtain agreement with the observations for the binary pulsar. The mean orbital 
radius for the binary pulsar is equal to the projected semi-major axis of the binary, 
{r)N = 7 X 10^° cm, so that {r)^ -C 14 kpc. Thus, G = Go within the experimental 
errors and agreement with the binary pulsar data for the periastron shift is obtained 
for Kn < 4.2 km. 

For a massless photon E = and we have 



d'u GM 



+ u = ?>^u\ (78) 

For the solar system using G = Gq within the experimental errors gives the light 
deflection: 

4GoMo , , 

= ^ (™) 

in agreement with GR. 



8 Galaxy Clusters and Lensing 

We can assess the existence of dark matter of galaxies and clusters of galaxies in 
two independent ways: from the dynamical behavior of test particles through the 
study of extended rotation curves of galaxies, and from the deflection and focusing 
of electromagnetic radiation, e.g., gravitational lensing of clusters of galaxies. The 
light deflection by gravitational fields is a relativistic effect, so the second approach 
provides a way to test the relativistic effects of gravitation at the extra-galactic level. 
It has been shown that for conformal tensor-scalar gravity theories the bending of 
light is either the same or can even be weaker than predicted by GR PH E21 • To rem- 
edy this problem, Bekenstein [21] has recently formulated a relativistic description 
of Milgrom's MOND model, including a time-like vector field as well as two scalar 
fields within a GR metric scenario. However, the time-like vector field violates local 
Lorentz invariance and requires preferred frames of reference. 
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The bending angle of a light ray as it passes near a massive system along an 
approximately straight path is given to lowest order in f^/c^ by 



J \a^\dz, (80) 



where _L denotes the perpendicular component to the ray's direction, and dz is the 
element of length along the ray and a denotes the acceleration. The best evidence 
in favor of dark matter lensing is the observed luminous arcs seen in the central 
regions of rich galaxy clusters jHS]- The cluster velocity dispersion predicted by the 
observed arcs is consistent within errors with the observed velocity dispersion of 
the cluster galaxies. This points to a consistency between the virial mass and the 
lensing mass, which favors the existence of dark matter. 
From ()78j) . we obtain the light deflection 

^^4<™^«^_ (81) 

c^R c^R 

where 

M = m(i + \^]. (82) 



M 

The value of M follows from for clusters as r ^ tq and 



G(r)-Goo = Go(l + ^^). (83) 

We choose for a cluster Mq = 3.6 x 10""^^ Mq and a cluster mass Md ~ 10^'* Mq, and 
obtain 




(84) 

We see that M ~ 7M and we can explain the increase in the light bending without 
exotic dark matter. 

From the formula Eq. (|56|) for r ^ Tq ~ 14 kpc we get 

«('•) = (85) 

We expect to obtain from this result a satisfactory description of lensing phenomena 
using Ea. pH) . 

The scaling with distance of the renormalized gravitational constant is seen to 
play an important role in describing consistently the solar system and the galaxy and 
cluster dynamics, without the postulate of exotic dark matter. In the next section, 
we outline a framework for an effective, running Newtonian gravitational constant 
G = G{k) and the skew field coupling constant 7c(fc), based on a renormalization 
group approach to gravity theory. 
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9 Renormalization Group Flow Gravity Theory 
and Running Constants 

We have postulated a renormalized Newtonian gravitational constant G = GqZ in 
our action Eq.Q, and in the phenomenology for galaxy dynamics 

(86) 

We shall now implement the running of the effective G and the -F^i^a field coupling 
constant 7c by using renormalization flow arguments in which the "classical" cou- 
pling constants G and 7c possess a scale dependent running behavior obtained by 
solving appropriate RG equations. This approach to quantum gravity and its astro- 
physical implications for the infrared (IR) low energy or large scale dependence of 
gravity has been applied to Einstein gravity [SI EH E3 tHH] with the intent to explain 
galaxy rotation curves without dominant dark matter. 

We shall postulate scale dependent effective actions rfc[(yf^,^] and rfc[y4^,^], corre- 
sponding to "coarse-grained" free energy functionals, which define an effective field 
theory valid at the mass scale k or length scale i = 1/k. Then the Tk are the bare 
actions obtained by integrating out all quantum fluctuations with momenta larger 
than an infrared cutoff kj^, or wavelengths smaller than i. The Tk are determined 
by functional differential equations for the RG flow, including the dimensionless 
coupling constants G{k) = k'^G{k),jcik), the cosmological constant X{k) = A(fc)//c^ 
and the mass fl = fi{k)/k. The dimensionless constant G{k)A{k) = G{k)X{k) is the 
quantity that is determined by measurement. The flow equations are described by 
a system of ordinary coupled differential equations: 

^^<5'(A;) = /5g(G',7c, A,/i), ^^7c(fc) = /^^-.(G', 7c, A, 

fc^A(A;) = /5A(G,7c,A,/x), fc^M^) = 7c, A, (87) 

To solve the problem, we are required to restrict the RG flow to a flnite-dimensional 
subspace P|, corresponding to a truncated theory space. In the truncated space 
only the coupling constants G and 7c, the cosmological constant A and the mass fl 
are considered. 

It is assumed that there exists an RG trajectory {G{k), 
7c(fc), A(A;), /i(fc)) and that this gives rise to running coupling constants G{k),-yc{k), 
a running cosmological constant A{k) and mass fi{k) with a cutoff k = k{x). The 
latter converts the scale dependences of G{k),'yc{k), A{k) and fi{k) to a spacetime 
position dependence jH[ EH EH] : 

G{x) = G{k = k{x)), 7c(x) = 7c(fc = k{x)), A = A{k = k{x)), 

fi = fi{k = k{x)). (88) 
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We shall follow the procedure of Reuter and Weyer [HI OH EHj and rewrite the 
action Sq in (j21) as 



and the skew field action as 

1 „ „„.,„ 1 



Sf = J d^xy^i-F^^.F^^'^P - -fx\x)A^,An. (90) 

We have required the skew field mass /i = ji^x) to be a function of the space- 
time coordinates, allowing for a mass renormalization. The matter action Sm is 
supplemented by an energy- momentum tensor 9^1,: 

' -ie„„, (91) 



which describes the energy-momentum associated with the G{x), 7c(a;), A{x) and 
fi{x). The new field equations that replace © are 

G^u + Ag^, = 8nGT^,. (92) 

The modified energy- momentum tensor T^^", which includes the contribution from 
B^'^, must satisfy the generalized Bianchi identities, which follow from the diffeo- 
morphism of the action: 

Vuf^" = 0. (93) 

These conservation laws lead to constraints on the G^'^. The fields G{x), 'ydx), A{x) 
and fi{x) are not varied in the effective actions Sq and Sp as they are treated as 
"background external fields" and they do not have corresponding field equations |Hl 

EH] . Thus, the actions do not contain kinetic energy terms corresponding to the 
fields (^(a;), 7c(a;), A{x) and /x(x). 

The flow equations are described by the running constants 

G{k) = G{k)/e, -f,{k), A{k) = Xk\ fi{k) = fi{k)k. (94) 

Every solution ((5(/c), 7c(/c), X{k),fi{k)) of the truncated flow equations is associated 
with the one-parameter family of action functionals 

^'^^^"^ = 16^^ / ^'^v^[^ - 2A(/c)], (95) 



TfclV] = / d'x^i-F.^^F^'^P - -^^\k)A^,An. (96) 

The effective average action formalism follows the RG flow from the bare action 
Ffc^oo = corresponding to the initial condition for the r^-flow equation, down to 
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k = and Tk=o = T corresponding to the effective action. We have to solve the 
effective equations of motion 



5T,[A 



0, (97) 



0. (9J 



In practice, we truncate the space of solutions and for the UV sector, simple 
local truncations are sufficient to describe the physical system, but in the IR limit 
/c — > nonlocal terms must be included in the truncation procedure. 

We shall presently just consider the running of G{k) and the coupling constant 
7c(fc). The RG flow equations are dominated by two fixed points and 7c*: Gaus- 
sian fixed points at = 7c* = and non-Gaussian ones with (5* > and 7c* > 0. 
The high-energy, short distance behavior of the quantum gravity theory is governed 
by the non-Gaussian fixed points, so that for — cx) all the RG trajectories run into 
these fixed points and lead to a quantum gravity theory by taking the ultra-violet 
(UV) cutoff along a trajectory running into the fixed points. This would lead to a 
non-perturbatively renormalizable quantum gravity theory |87 | l88 t l89 t I41j . The 
conjectured existence of such a non-perturbatively renormalizable gravity theory 
still has to be demonstrated. 

Both the couplings G{k) and 7c(fc) are asymptotically free coupling constants 
as in quantum chromodynamics (QCD), so that they vanish for k oo. We see 
that if we interpret i = l/fc as a distance scale, then both G{k) and 7c(fc) increase 
all the way from the UV energy scale to the infrared (IR) energy scale as — >■ 0. 
This corresponds to a quantum gravity anti- screening, playing the analogous role to 
the anti-screening in Yang-Mills QCD. The RG trajectory ending at the fixed point 
G{k = 0) = can be described by the scaling law 

G(fc) = Go(l + ^). (99) 

This interpolates between the constant Newtonian value G{k\ab) = Gq at large k 
corresponding to short distances and the increasingly larger values of G(k) at larger 
distances for small k. The difference between G at k = and at a typical laboratory 
scale, k ~ l/(meter) is negligible. Any UV renormalization effects are negligible for 
laboratory, galaxy and cosmological scales. 

The gravitational coupling constant will display a running power law behavior: 

G{k) = ^. (100) 

For the spherically symmetric static solution, the running of G{r) is described by 
Eq. (j3H|) . so that as r increases to galactic distances G{r) increases until it reaches 
the constant value 



Goo — Go 




(101) 
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A more detailed analysis of RG flow equations in terms of our effective actions 
will be presented elsewhere. 

We note that if we had considered an effective RG framework based on the 
Einstein-Hilbert action Q without including a coupling to the skew field F^j^x, then 
the weak gravitational field acceleration law would take the effective "Newtonian 
form" (j57jl . Then, the RG flow effective action formalism is required to predict that 
G{r) should run in such a way that it produces flat rotation curves for galaxies 
that flt all the galaxy data [HI]. Whether this is possible is a matter of conjecture. 
In contrast, the running of the effective constant G = G{r) in MSTG theory is 
determined to a large extent by the fleld equations of the theory for weak flelds and 
leads to good agreement with the available galaxy rotation curve data. 



10 Cosmology Without Cold Dark Matter 



We shall assume that the universe is isotropic and homogeneous at large scales and 
adopt the FLRW line element 



ds^ = de -R\t) 



1 — kr"^ 



The Friedmann equations take the form 



47rG 



InGpMjt) A 
3 3' 



A 



[pM{t)+3pM{t)]R{t) + -R{t), 



(102) 



(103) 



(104) 



where H{t) = R{t)/R{t). 

For a homogeneous and isotropic universe the skew symmetric flelds A^,^ and F^j^a 
are zero when averaged over large distance scales, since there can be no preferred 
direction in the maximally symmetric FLRW spacetime. In Eqs. ()103|) and ()104|) . 
we have assumed that G = GqZ and the renormalized value of G is obtained from 
the running of the effective coupling G = G{t) with time. We can write p()3j) for a 
spatially flat universe = in the form: 



where 



SirGpMit) 



A 



(105) 



(106) 



We shall impose physical restrictions on the form of the running of G{t). We observe 
that from calculations of the production of helium and deuterium abundances at 
the time of BBN ^42j, it can be shown that fi_B ~ 0.02 — 0.04. Moreover, from 
the WMAP data the fractions of baryon and photon densities at the surface of last 
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scattering show that Qb ~ 0.04 Therefore, we must impose the conditions 
that at the time of nucleosynthesis at tBBN ~ lO^'^yrs (at red shift z ~ 10^), 
G{tBBN) ~ Gq and at the surface of last scattering (decoupling) at tsis ~ lO^yrs 
[z ~ 10^), we require that G{tsLs) ~ Gq. Moreover, since the time of the surface of 
last scattering we additionally require that G{t) grows until the present time to the 
value G{t 

now) ~ QGq. Given this running of G, it follows that for t > tgi^s'. 

VLm ~ Q^B ~ 0.24, (107) 

and we need only assume a dominant baryon density. Given this scenario, we do 
not require a dominant cold dark matter in the matter dominated era. There can 
be small contributions to VLm due to massive neutrinos. 

A bound on the possible value of a changing G with time is obtained from recent 
spacecraft measurements |43j : 

G 

- = (4±9) X lO^^^yr"^ (108) 

Gr 

In order to apply this bound to determine the amount of variation of G(t), we need 
to know the functional form of G{t) as it runs with time. This must be determined 
by solving the RG flow equations for our cosmological model. If the running of 
G with time is such that G/G ~ for a flat behavior of G{t) ~ Gq up till the 
surface of last scattering, followed after the surface of last scattering by a rapid rise 
to G{t) ~ 66*0 with a subsequent value GjG ~ 0, then the bound p08p will not 
rule out the required renormalized value G(tnow) ~ QGq. In Fig. 6, we display a 
conjectured effective G{t) as a function of t in years that is consistent with (jl08p 
and a baryon dominated universe without exotic cold dark matter. 
The value oiG{t 

now) ~ SGq required to obtain results consistent with the WMAP 
data 7j corresponds to cosmological distance scales. For the solar system with a 
distance scale of order 1-40 AU, the local value of G'(tnow) ~ Gq to within the 
experimental errors on the measurements of Go- 
lf we take into account the running of the cosmological constant A (A;) in the RG 
quantum gravity flow framework, then we can have A(t) run with t in ()103p and 
fjl04j) and this would correspond to a quintessence scenario jUj. It is possible that 
the RG flow trajectories that lead to a large distance classical scenario can solve 
the cosmological constant problem, for these trajectories imply a small cosmological 
constant [311 EE] • 

A comparison of our baryon dominated matter era cosmological model with the 
acoustical peak data obtained from the WMAP observations [7] has to be performed 
to see whether a satisfactory fit to the data can be obtained. The concordance 
KG DM model agrees well with the data and the present model must produce fits 
to the CMB data that succeeds as well. Moreover, an investigation of the growth of 
large scale galaxies and clusters from an initially smooth background of fluctuations 
must also be performed. Since the effective G in the cosmological model is required 
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to increase after the surface of last scattering to a value G{t 

now) ~ 6Go, then we 

can expect that the gravitational potential well will become deeper and allow for 
clumping of the baryon matter to form galaxies without cold dark matter. The 
results of this investigation will be reported elsewhere. 

11 Conclusions 

We have developed a gravity theory consisting of a metric-skew-tensor action that 
leads to a modified Newtonian acceleration law that is fitted to galaxy rotation 
curves. There is a large enough sample of galaxy data which fits our predicted 
MSTG acceleration law to warrant taking seriously the proposal that the gravity 
theory can explain the fiat rotational velocity curves of galaxies without exotic 
dark matter. It is interesting to note that we can fit the rotational velocity data 
of galaxies in the distance range 0.02 kpc < r < 70kpc and in the mass range 
10^ Mq < M < 10^^ Mq without exotic dark matter halos We have also 

provided a model for the behavior of galaxy rotation curves inside the luminous 
core of galaxies that predicts well the observed flatness of the rotation curves as the 
distance from the core increases. The lensing of clusters can also be explained by 
the theory without exotic dark matter in cluster halos. 

Our RG flow effective action description of MSTG quantum gravity allows for a 
running of the effective G with distance. The RG flow framework for the theory is 
characterized by special RG trajectories. On the RG trajectory, we identify a regime 
of distance scales where solar system gravitational physics is well described by GR, 
which is contained in MSTG as an approximate solution to the field equations. We 
are able to obtain agreement with the observations in the solar system, terrestrial 
gravitational experiments and the binary pulsar PSR 1913+16. Strong infrared 
renormalization effects become visible at the scale of galaxies and the modified 
Newtonian potential replaces exotic dark matter as an explanation of fiat rotation 
curves. Thus, gravity becomes a "confining force" that has significant predictions 
for astrophysics and cosmology. 

We have demonstrated that the RG flow running of G and MSTG cosmology 
can lead to a description of the universe that does not require dominant, exotic 
dark matter. Dark energy is described by an effective time dependent cosmological 
constant. A detailed investigation of the MSTG cosmological scenario must be 
performed to establish that it can describe the large scale structure of the universe, 
account for galaxy formation and big bang nucleosynthesis and be consistent with 
the WMAP data. The RG flow trajectories in MSTG in the IR large classical limit 
can lead to a possible solution to the cosmological constant problem, for they imply 
a small cosmological constant [SI EH] • 

The skew fields F^^^cr have a natural interpretation in terms of string theory ^T] , 
so that a possible string theory interpretation of the astrophysical predictions can 
be investigated. 
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An interesting approach is to solve the MSTG field equations for gravitational 
collapse to see whether there are any stable massive skew field solutions that cor- 
respond to conventional Schwarzschild black holes with event horizons. If the con- 
jectured asymptotically safe renormalizability holds for MSTG quantum gravity, 
then we might expect that G{s) and 7c(s) vanish as the distance s = {x ■ xY^'^ — > 0, 
which could lead to a gravitationally-free Minkowski spacetime without singularities 
at s = 0. 
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Table 1. Values of the total galaxy mass M used to fit rotational velocity data. 
Also shown are the mass-to- light-ratios M/ L with L obtained from ref . ^\ . 
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Fig. 1 - Fits to low-surface-brightness and high-surface-brightness spiral galaxy 
data. The black curve is the rotational velocity v versus r obtained from the 
modified Newtonian acceleration, while the green curve shows the Newtonian 
rotational velocity v versus r. The data are shown as red circles. 
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Fig. 2 - Fits to the data of two dwarf galaxies Fornax and Draco. The simple 
relation V ~ ^/2a is assumed between the velocity dispersion a and the rotational 
velocity v. The black curve is the rotational velocity v versus r obtained from the 
modified Newtonian acceleration, while the green curve shows the Newtonian 
rotational velocity v versus r. The data are shown as red circles and the errors 
(not shown) are large and for Draco the Newtonian fit cannot be distinguished 
from the MSTG fit to the data within the errors. 
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Fig. 3 - Fits to the elliptical galaxy NGC 3379 and the globular cluster u 
Centauri. The black curve is the rotational velocity v versus r obtained from the 

modified MSTG acceleration, while the green curve is the Newtonian-Kepler 
velocity curve. For u Centauri the black and green curves cannot be distinguished 

from one another. 



30 



Fig. 4 - Fits to the rotation curves for the four galaxies NGC 1560, NGC 2903, 
NGC 4565 and NGC 5055. The black curve is the rotational velocity v versus r 
obtained from the modified MSTG acceleration formula (j65|) . while the green curve 
is the Newtonian-Kepler velocity curve for the core-Newtonian acceleration 
formula (jU^ . The luminous mass used to fit the data for NGC 1560 is 
M = 1.51 X 10^° Mq, while the core radius is Vc = 1.5 x 10^^ cm = 4.85kpc. For 
NGC 2903 the luminous mass used to fit the data is M = 9.95 x 10^° Mq, while 

the core radius is Vc = 2.75 kpc. For NGC 4565 the luminous mass is 
M = 1.81 X 10" Mq, while the core radius is = 3.11 kpc. For NGC 5055 the 
luminous mass is M = 8.55 x 10^° Mq and the core radius is = 2.04 kpc. 
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Fig. 5 - 3- dimensional plot of v versus the range of distance 0.1 kpc < r < lOkpc 
and the range of galaxy mass 5 x 10® Mq < M < 2.5 x 10^^ Mq. The red surface 
shows the Newtonian values of the rotational velocity v, while the dark surface 

displays the MSTG prediction for v. 
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Fig. 6 - Plot of conjectured running of the effective gravitational constant G{t) 

with time t. 
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